Twist-three analysis of photon electroproduction off pion. 
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t:;^- ' Abstract 

cn 

We study twist-three effects in spin, charge, and azimuthal asymmetries in deeply virtual 

O Compton scattering on a spin-zero target. Contributions which are power suppressed in 1/Q 
O ' 

' generate a new azimuthal angle dependence of the cross section which is not present in the leading 

Q^i twist results. On the other hand the leading twist terms are not modified by the twist three 

^ contributions. They may get corrected at twist four level, however. In the Wandzura-Wilczek 

'-^ approximation these new terms in the Fourier expansion with respect to the azimuthal angle are 

• ^ ', entirely determined by the twist-two skewed parton distributions. We also discuss more general 

>< ; ..u. Uke the gene.al fonn of the ^gula. dependence of the difletential cto. sectton, valtdtty 

■ ■ ■ ' of factorization at twist-three level, and a relation of skewed parton distributions to spectral 
functions. 
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1 Introduction. 



Hard electroproduction of a photon, i.e. deeply virtual Compton scattering (DVCS) |T], ||, 
measures non-forward matrix elements {P2\Oi\Pi) of quark and/or gluon non-local composite op- 
erators Oi- The former possesses a more diverse structure than processes with forward kinematics 
e.g. arising in conventional inclusive reactions. Fourier transforms of {P2\Oi\Pi) define a new type 
of hadron characteristics — generalized parton distributions (GPDs). Real photon electroproduc- 
tion on nucleon targets has already been studied by DESY experiments: for moderate Bjorken 
varaible x-q by HERMES ||^ and for small by ZEUS and HI collaborations 0. Since these new 
functions are poorly known, our primary goal is to constrain diverse models from experimental 
data. To extract GPDs from future experiments |^ one has to address the problem of appropriate 
observables. The four-fold differential cross section for e{k)h{Pi) e{k')h{P2)'j{q2) 



da a^XBV _^4:M'^x'^\ 



dxBdyd\A^\d(l) 8 7rQ2l q2 



depends on the Bjorken variable Xb = —q\/{2Pi ■ qi) (with qi = k — k' and q\ = — Q^), the 
momentum transfer = (Pg — -Pi)^, the fraction of the lepton energy loss y = Pi ■ qi/Pi ■ k 
and the azimuthal angle (p between lepton and hadron scattering planes. The amplitude T is 
a sum of the virtual Compton scattering (VCS) amplitude, Tycs: and the Bethe-Heitler (BH) 
amplitude, Tbh, since they have the same initial and final states. Fortunately, the interference 
term, X = Tycs'^bh'^'^vcs'^bh, in I'^P may be isolated by measuring various asymmetries like the 
lepton charge asymmetry, the hadron/lepton spin asymmetry and the azimuthal angle asymmetry 
[0, 1^, 1^. Since X is linear in the DVCS amphtude, one can extract its real or imaginary part 
from an observable defined in an appropriate manner. Thus, one directly accesses a specific linear 
combination of GPDs convoluted with the real or imaginary part of a hard scattering amplitude. 

The theoretical description of DVCS amplitudes is simplified for small values of xb- This region 
is relevant for the HERA collider experiments. For small xb one naively assumes that GPDs are 
essentially determined by the forward parton densities^ multiplied by partonic form factors. The 
momentum transfer dependence of GPDs may be parametrized by a hadron form factor or for 
larger values it can be computed within QCD |jlO|- For larger values of xb the 'exclusive' domain 
of GPDs becomes important. Then their functional form cannot be approximated by conventional 
densities. Rather one has to rely on models or plausible parametrizations Such predictions can 
be unstable under QCD radiative- and/or higher twists corrections. Therefore both issues deserve 



a detailed investigation. Although the evolution effects are quite moderate (of order 10 — 15% [|TT|), 
the next-to-leading order corrections in the coefficient functions have a significant impact on the 



^The only free parameter is the renormalization group scale where this identification is done. 
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handbag approximation and can be of order 50% in the valence region |T^. On the other hand the 
exploration of power suppressed contributions has just begun. Recent studies explored twist-three 
GPDs in the DVCS amplitudes ||13|-|17] and relations of twist-three GPDs to already discussed 



leading twist functions [|T^]-||T3 and interaction dependent antiquark-gluon-quark correlations [|T^. 
A natural next step is to study power suppressed contributions in the coefficients of the Fourier 
expansion of the differential cross section. For the sake of simplicity and clarity we consider a pion 
target, h = n. While barely a target for an experiment, (pseudo) scalar particles are instructive 
for theoretical reasons. They avoid difficulties accompanying particles with non-zero spin. Our 
goal is thus to gain insights into gross features of twist-three effects in physical observables which 
are expected to hold true also for DVCS on higher spin targets. Presently, we discuss issues of the 
azimuthal angle dependence of the differential cross section (|l|) with contributions from twist-three 
GPDs and qualitative estimates of higher twist effects. 

The paper is organized as follows. In Section |^ we calculate the azimuthal angle dependence 
of the differential cross section in terms of the hadronic amplitudes, which appear in the DVCS 
tensor for a spin-zero target. We derive constraints between the angular moments of the amplitude 
squared. In Section ^ we evaluate the twist-three predictions based on perturbative leading order 
calculations. In Section ^ we give a consistent definition of spectral functions, so-called double 
distributions, and derive their relationship to GPDs. After giving a few numerical estimates in 
Section ^ we summarize. 



2 Angular dependence of the cross section. 

As we observed before, the total electroproduction amplitude of a real photon consists of two 
terms, Tycs, "we are interested in, and the contaminating Tbh- Tbh is purely real and arises 
from a contraction of the leptonic tensor, L^^ = u{k', A') [j^i^ + 7i.(|^' +/^)^^'^^ u{k, A), 

with a hadronic electromagnetic current, which for the pion target reads = P^F(A^). It is 
parametrized by the electromagnetic form factor F(A^) described by a simple monopole form 
(1 — A^/my) ^ with my ~ 0.46 GeV^. Similarly, the DVCS amplitude is expressed by Tycs = 
■^^T^ye*^{q2)u{k')^yu{k) , where the sign depends on the beam charge, ± for e^. The hadronic 
tensor is given by 

T^,(g,P, A) = I J dxe'^-SP2\T3,{x/2)3,{-x/2)\P^). (2) 

It depends on the momenta q = (gi + g2)/2, P = P1 + P2 and A = P2 — P1 which can be reexpressed 
in terms of the Lorentz invariants ^ = Q^/P ■ q, f] = A ■ q/ P ■ q, = —q^ = —j{qi + q2Y 
and A^. For the DVCS kinematics, the first two scaling variables are related to each other by 
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1] = (^1 — ^ —C,. The 'experimental' variables referred to in the introduction are expressed 
by the present ones to twist-four accuracy via, see e.g. ~ ^ ~ 3;b/(2 — xb). There are 

five independent kinematical structures in Eq. for a general two photon process on a (pseudo) 
scalar target: 

■ q 2F ■ q 

+ 2p-q + IP 

Here current conservation is ensured by means of the projector V^^ = g^u — qi^q2u/<li ■ 12- The 
transverse component of the momentum transfer is = A^ — rjP^. We scaled the kinematical 
factors in such a way, that all five dimensionless scalar amplitudes change only logarithmically if 
quark masses are set to zero. Contracting the hadronic tensor with the leptonic current we can 
conclude that the amplitudes Fi, F2 and F5 are the leading contributions in Q, while F3 and F4 
are suppressed by the powers 1/Q and respectively. Note also that the Bose symmetry 

requires Fi, F2, F4, F5 to be even and F3 to be odd in 77. For the case at hand the outgoing 
photon is real. Then the on-shell condition ^2 = forces V^^^a to be proportional to ^2/^- When 
this term is contracted with the leptonic part its contribution vanishes. Only three contributions 
are left then 

:i-a(A^- Ag 



T^.(g, P, A)| = -V,.g.rVru T, + ^4—- ^T, 



V P P V / \ VP A-^V 

+ 2P.g 1,^^ + 4M^^V 2M^ 

The new form factors are parametrized in such a way that gluonic transversity only contributes 
to T3 at 0{as). We also dropped T^jQ? contributions in the first Lorentz structure. The new 
amplitudes are parametrized by Compton form factors, which are related to the old ones by the 
following set of equations 

\ / |r7=-5 V / \'q=-i 

We again neglected twist four terms. Due to current conservation, the substitution P^ H 
may be performed. It is obvious then that the second term on the r.h.s. of Eq. (^) is suppressed 
by at least one power in 1/Q w.r.t. the leading terms which contain Ti and T3. Once more we 
emphasize that is absent at the Born level. 

Now we are in a position to discuss the azimuthal angle dependence of the DVCS and BH 
amplitudes squared as well as of the interference term. To ensure that the Fourier series con- 
tains no terms which are artificially generated by kinematical subtleties, we choose a frame 
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rotated w.r.t. the laboratory oneQ. In our frame the virtual photon has no transverse com- 
ponents and a negative z component q\ = {q1,0, — \qf\). The x component of the incoming 
positron is positive k = {E, E sm6e,0, E cosOe). Other momenta read Pi = (M, 0,0,0) and 
P2 = (-E2, 1-P2I cos0sin6'i:/, IP2I sin0sin6'H, I-P2I cos6'/f). is the azimuthal angle between the 
lepton and hadron scattering planes. 

The calculation of the squared amplitudes is straightforward and yields (the electric charge is 
set equal to one): 

• Bethe-Heitler squared term. 

The amplitude squared for a spin- zero target does not depend on the lepton polarization. It reads 
,2 8F2r4M2-A2/ Q^ + A^ \ {{q - k) ■ Pf + {k ■ Pf \ 

•^BH — —TTT { TTTT, 1 + TTT^t; 1 . on / . ^, 7T + 



A2 [ 2A2 \^ 2(2A:- A- A2)(Q2 + 2A;- A)y • A - A2)(Q2 + 2A; ■ A) J 

(5) 

The denominators of both lepton propagators {2k ■ A — A^) and {2k ■ A + Q^) depend on the 
azimuthal angle 0. Moreover the Taylor expansion of {Q"^ + 2A; ■ A) = {k — ^2)^ = + 
0{1/ Q)) induces a pole for y ^ 1. We introduce as a new notation Vi and V2 for the dimensionless 
lepton propagators: 

Q^'Vi = {k- q^f = Q^ + 2k-A, Q^V2 = {k- = -2k ■ A + A\ (6) 
Their dependence is contained in 



where e = 2x-bM/ Q. The kinematical factor 
K = 



'V 4 ;v / r 4(1 -xb) q2 

is l/Q power suppressed. It vanishes at the kinematical boundary A^ = A^j^^ determined by 



-^-^--^ 4xB(l-XB)+e2 -1-xb + xbM2/Q2|1 + ^1^^/2 j|- (9) 

Let us take a closer look at equation (||). In the second numerator, {{q — k) ■ Py + (A; ■ P)^ = 
{q ■ P — k ■ A — 2k ■ PiY + {k ■ A + 2k ■ Pi)^, only k ■ A depends on 0. Therefore equation (|^) 
contains terms proportional to cos° 0, cos^ and cos^ 0. In other words (|^) may be written as 

•^This reference frame is related to the centre-of-mass system in Ref. Q by a boost of the hadron in the z- 
direction. 
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where the expansion coefficients are given by 

C = {{2~y? + y\l + |4x|^ + 4(1 - xb) + (4xb + 

+32x1^2^ + 2e' {4(1 - y){3 + 2e') + y\2 - e^)} - 4x|(2 - y)2(2 + e^)|!, 
cfH = -8(2 - y) |2xB + - 4x|^| , 

cf = 32x|^. (11) 

It should be noted, that all coefficients are bounded in the whole kinematical region, since ^ is 
always multiplied by a;| and A^j^ oc M^x|. The K term in the Fourier coefficients Cq, Ci and 
C2 suppresses the angular moments by / rf^cos (m0)|7BHp { ^ ~ %^)} ^ • '^'^ prevent 
induced azimuthal angle dependence due to the expansion of lepton propagator in l/Q, angular 
moments with an aditional weight of the two propagators V1V2 can be calculated. Then the 
averaged BH squared term reads for large 

f '4v,(Mr.uP ^ {8(2 - 2y + ,^)(1 - .e) (l - + O (i,) } . (12, 

• Bethe-Heitler-DVCS interference term. 
The interference term can be treated in the same way. 

^ = - A2ri.^2^ l^^o + E \cl COS (m0) + Xsi sin (m^)] + ^iT^cf cos (30)| (13) 



m=l 



As we see the Fourier sums terminate with cos (30) and sin (20) terms for unpolarized and polarized 
scattering, respectively. The information about GPDs is contained in the dimensionless coefficient 
functions and s^. They are linear combinations of the amplitudes Tj introduced in (1). The 
exact results are rather lengthy. Therefore we neglect terms which are power suppressed w.r.t. 
the leading contribution: 

= -8xb(2 - y) |(2 - xb)(1 - y) - (l - xb)(2 - y f (^1 - | ReTi, 

4 = -8xB{2-2y + y^)ReT,, 

d = ^ReTg, 

' 2 - xb 

= 8xBy(2-y) ImTi, 

4 = 16ylm |xbTi - - ^Tgl . (14) 
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As we can see, the interference term falls off like Ij \J —IS^Q} for large Q. Its average over 
is suppressed by Note that the coefficient cf arises at twist-two level due to the gluon 

transversity in T3. 

• DVCS squared term. 



Finally, the DVCS squared term reads 



TDvcsr = \ cr^' + K cr^^ cos (0) + a .^^^ sin (0) + ^i^^c^^^ cos (20) \ . (15) 
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^DVCS /' I \ „DVCS 



2 



-2^DVCS 



M2' 



The Fourier expansion terminates with cos (20) and has only a sin (0) term. The coefficients are 
quadratic in the Compton form factors Tj and to leading order in 1/ Q they are given by 

= 2(2 -2y + y-) {xbT^T* + ^^j^^^^^^,§^T,A , (16) 



'1 




2;b(2 - x^Y 
A2 - A2 . 

mm rri 

2-xb a;BM2 ^ 1"^"^' 2 - xb 



cr^^ = 8(2 - y)M ( - V""^^ :tr TA \ xbT, - T, - l^^^i^T, 

^DVCS 

z — XB 

nvcs _ o„,T_ J ( ^ 1 - xb - A^i„ \ / ^ ^ 1 - xb A^ - A^j^^, 



xbTi - T2 - — 



2-xb a;BM2 A 2 - xb 

The Compton form factor T3 contributes to c^^'"^ and c^^*"^ at leading twist. In the remaining 
two coefficients the interference of Ti with T3 (as well as in cf and sf) requires a twist-three 
analysis at next-to-leading order. Note that all angular coefficients are determined by geometrical 
(i.e. dimension minus spin) twist-two (for cf , cf , sf , cf , Cq^*"^ and c^^^^) or by geometrical twist- 
three contributions. However, cf and c^^*"^ may be stronger contaminated by power suppressed 
contributions than the other coefficients, since gluon transversity is perturbatively suppressed. 

It is important to note, that some of the coefficients are related to each other. We recall that 
there are only three Compton form factors in equation (^. Therefore, the four coefficients cf are 
reduced to three independent ones. To leading order approximation in 1/ Q this relation looks 
like: 




cf = 2-2yly^ {^^ " "^^^^ " ^) " " ^b)(2 - vY 1^1 ' ^ ] \ + O (i/Q^) . (17) 

Measuring charge asymmetries yields information about the real part of all three amplitudes and 
the imaginary part of two of them (see discussion below). If we assume for a moment that the T3 
contribution is negligibly small, we find the following relations 

, 2 



^ 2-2y + y^ (4) {s^ 

° 32XB (2-2y + y2)2 + ^2(2_^)2 ^+ I /-^ 
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DVCS _ ^ ) I I + o d/O'^) 

„DVCS -'- f ^1*^ '^■^2 I I /n /'l 



^fcf 

16xB 1(2 2-2y + y2 



+ 1/Q- 



These constraints can be used to test Q scaling of the Compton form factors. They can be 
generahzed to the case of non-zero T3. The real part of T3 enters in cf , while the imaginary part 
of T3 can be extracted from c^^*"^ (assuming we know Ti). In this case we obtain the same set of 
equations (|18]), but the coefficients cj^^*"^, cf^*"^ and sf^*"^ then depend on all angular moments 
of the interference term as well as on c^^*"^. 

Finally, let us give a few remarks concerning the actual measurements of the angular coeffi- 
cients. Of course, experiments with pions are improbable at present but the remarks also apply 
to the analogous case of lepton-nucleon scattering. For optimal experiments one needs both po- 
larized electron and positron beams. Since the amplitudes squared are charge-even, while the 
interference term is charge-odd, this allows to access all moments of the interference term. The 
cos (0) and sin (0) dependence should dominate away from the kinematical boundary, because 
other contributions are either kinematically or perturbatively suppressed. After subtraction of 
the BH squared contribution measurements of the charge-even part of the cross section for unpo- 
larized lepton beams allow to access all azimuthal angular moments of the DVCS squared term. 
Here the constant term is dominant, because K suppresses the other terms kinematically and 
because of the expected smallness of transversity. With a polarized beam the sin(0) term can be 
directly accessed, because the BH process is independent of the beam polarization. Following the 
outlined schedule one should, thus, be able to test the three equalities (pISl). 

Let us now comment on experiments where a polarized lepton beam is only available for 
one charge. Then single spin asymmetries may be used to yield information about the angular 
coefficients by weighting the differential cross section with ViP2- Unfortunately the DVCS squared 
term will inevitably contribute to the azimuthal asymmetries then. Although its contribution is 
suppressed by yA^/ [compare Eq. (|1^) with Eq. (p^S])] compared to the interference term it 
is also affected by the ratios ReTj/F. For the unpolarized case the following observations are 
relevant. Near the kinematical boundary = A^j^ all terms of the amplitude squared tend 
to be of the same order starting with Away from A^ = A^j^ one still has a relative 

kinematical suppression of ITdvcsP by y^'^/Q^. However, yA'^/Q^ is now multiplied by terms 
containing Im Tj/ F. Of course, one has to subtract the BH processes to get the desired information. 
Fortunately, for the cos(0) and cos(20) moments the BH term enters with the same power 1/Q 
as the interference term [compare Eq. ( |10|) with Eq. (|13])]. The latter again includes the ratios 
ReTi/F. 
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3 Twist-three GPDs. 



In this Section we derive the exphcit expression for the Fourier coefficients in terms of GPDs 
from the hadronic tensor, which is known to twist-three accuracy in the Born approximation. We 
will discuss different aspects of the result. Since factorization has not yet been proven for 1/Q 
suppressed contributions it might happen (in contrast to the leading twist situation 
that the results are plagued by singularities. Universality of higher twist distributions could not 
be guaranteed then. A keyrole in factorization plays the property that non-perturbative functions 
do have specific analytical properties. For instance, in the case of the pion transition form factor 
measured in 7*7* ^ vr*^, it requires that the meson distribution should vanish at the end-points 



and this can indeed be proven p2[. In the case of DVCS an analogous requirement is the continuity 



of skewed parton distributions at |x| = rj. The tree level analysis has been recently done |13|-[17|. 
For transverse polarization of the initial photon a divergent amplitude is encountered. However, 
this divergence does not show up in the cross section because of the final state 7-quantum e*-q2 = 
WTl. Note that this cancelation of singularities has not yet been proven for twist contributions 
higher than twist-three. 

For any target the hadronic tensor reads (e^^^s = x) in the considered order 



^ (a P A) = -V a V j-f-p p-p j^-p p -p ) Yl£ P ie V 

. q r- f F ■ q F ■ q 



Vip and Alp are given as matrix elements of vector and axial- vector light-ray operators, respectively. 
We should note that the gauge invariance beyond the twist-three accuracy has been restored by 
hands. Furthermore, in our approximation the vector V2p depends on the other ones by the 
relation: 

V2p = (Vip- ■ Vi + '/-^Ai^ + twist-four. (20) 

2 F ■ q 2 F ■ q 

The Lorentz structure of the Ai a- term in ( PP] ) when combined with the last term in (|T^) changes 
the gauge invariant V^^e^rqpVru projector into {g^^ - F^ q2a/P ■ ^2) ^arqp {Qut - Py Qir/P ■ Qi) [HI- 
The result ( p!9| - |20|) contains relations that are well known from deep-inelastic scattering. Obviously, 
in the forward case, where Vip oc PpFi and ^ = xb, we rediscover the Callan-Gross relation. Besides 
the generalization of this relation, there are other ones at twist-three level. They will be discussed 
below. In our case of a scalar target, the matrix element of the axial- vector operator has a non- 
zero form factor at twist-three level, while the vector operator already has a non-zero form factor 
at leading twist level: 

Vi p = Ppn + ( - T^Fp)n''^'' + twist-four, Ai p = H'^^^ (21) 
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The (twist-two and three) form factors Ti and Ti are given as convolutions in momentum fraction 
X of perturbatively calculable hard scattering parts with GPDs ((g) = / dx): 




(e,r^,A2|Q^)= E (e,a:|QV/^')® ~ (x, r/, A^j/i^). (22) 




n 

The summation runs over different parton species and /i^ is the factorization scale. To simplify 
notations we have kept only scaling variables as an argument of the functions. For a Qrcharged 
quark the leading order coefficient function reads ^ C^^'' (^,x) = Q}{1 — x/^ — it)^^ ^ {x —x) 
with — (+) standing for parity even (odd). 



The complete twist decomposition of the vector and axial-vector operator is given in ||15[ and 
our parameterization ([21| ) reads for the twist-two GPDs 

{h,h) {x,v) = J ^e-^+^'(Ar'^0+(«:,-«:)|Pi), ^'^0+(k, -«:) = V5(-Kn)7+ (1, 75) ^(«:^).(23) 

We omitted gauge links here. Note that these functions can be related to other fundamental 
non-perturbative amplitudes, namely hadronic wave functions, via overlap integrals [|18|. The 
twist-three functions contain two pieces, the so-called Wandzura-Wilczek (WW) piece expressed in 
terms of the twist-two function and a correlation function, G, of antiquark-gluon-quark operators 



15 



H'-'i-^v) = /; - 1 I I -| I Hiy,,) + l H{x,v) - G+(x,,), 





1 X 


y 


-1 \r]\ 


\V 


V 




/ X 


y 


-1 \ri\ 


\V 


V 



Here W±{x, y) = ^ {W{x, y) ± W{—x, —y)} and the W kernel reads W{x, y) = 6'5'i(l + x,x — y) 
with |]T^ Oii{x,y) = Jq da6{xa + ya) = [0{x) — 6{y)]/[x — y]. Note that the first moment with 
respect to x vanishes. Moreover, the time reversal invariance together with hermiticity tells us 
that the twist-three GPDs are real valued functions with the symmetry properties if*™^(— r^) = 
— if*™^(?7), if*™^ (—■?]) = i7'™^(?7). The only new dynamical information is contained in the 
antiquark-gluon-quark GPDs 



(25) 

We used the convention W" = ^iy(±-,±-). Due to this second derivative the 

antiquark-gluon-quark operators do not contribute to the second moment of twist-three GPDs. 
Since the operators defined in the parity even and odd sectors can be expresses through each other. 
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only two independent functions {P2\S'^{n,UK,—K)\Pi) = P+A^ / dx exp{—iKxP+)S^{x,u,ri) ap- 
pear, where the operators are iS^(ki, K2, ^^s) = ig^{i^zfT) l+G+p{i^2n) ± i7+75G'+p(K2n) ip^Kin). 
Note that the twist-three form-factors possess discontinuities |jl6|, [l^ at the points x = \ri\, for 
instance: 

H^'^^T] + 0,v)- H'-^Hv - 0, v) (26) 

where we assumed that the GPDs vanish at the boundary x = ±1. The appearence of disconti- 
nuities is a general artefact of the procedure to separate twist-two and -three contributions and is 
not solely due to the WW- approximation. 

There is a number of relations between the amplitudes Fi in at the Born level. The first 
one is a generalization of the twist-two Callan-Gross relation: 

F,i^,v) = ^F2i^,v)= E f dx C\-\x,i)H,{x,ii). (27) 

^ i=u,d. ^ 



This equation can be easily deconvoluted by J2i [Hi{x, rf) — Hi{—x, r])] = ImFi(x, ri)/iT. Moreover, 
by means of antisymmetry we can continue Fi{x,ri) to negative values of x. Thus, we obtain a 
dispersion relation, vrRe Fi(,^, r/) = PV j\ 7^ImFi(x, r/), which tests the dominance of perturba- 
tive leading order predictions. In the WW approximation for twist-three functions two further 
relations hold, which express the amplitudes F3 and F4 in terms of twist-two GPDs: 

p - ^ Q'i ( sig^(^) i + y-i^ +{r] -v)] ( — — 

' i=ud * \\?7 + x ^ ~ X — iO / 1] dx 

2 




■r]{^ — X — iO) 

m.) - E «£-{ (^-f^ - (. - -.)) (I ~ i 

-{x^-x)^Hi{x,7]). (28) 

As we can see, in the kinematics of a two-photon process the WW approximation exist, although, 
the ln(^ ± 7]) terms could provide a numerical enhancement. In the case of an outgoing on-shell 
photon, i.e. ^ = — r/ + (!?(A^/Q^), one naively expects a logarithmic enhancement (stemming from 
(p6D ) due to ln(^ + ^7) ln(A^/ Q^), which would imply difficulties with factorization. However, it 
has been shown that such contributions cancel in physical amplitudes [ll7| , 0, and thus legitimize 



the WW approximation. It is interesting to note that the same situation appears also in the 
antiquark-gluon-quark sector. 
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The DVCS (77 = —^) amplitudes in the tensor decomposition (|[) are given by 

Ti =n, T2 = 2^n + |7^tw3 _ ^tw3| ^ (29) 

and the form factor T3 is determined at twist-two level by the gluon transversity. All Compton 
form factors also contain twist-four contributions, which have not been computed yet. Obviously, 
the generalized Callan-Gross relation can not be tested in the DVCS process, since it is modified 
by twist-three contributions, see (^91). The modification is given as a difference of 7i*™^ and 7i*™^ 
and, thus, ensures the cancellation of the ln(^+?7) term for the WW and the antiquark-gluon-quark 
contributions: 



where Si{x,^) = j\ du{l + u) Sf{—x, —u, — ^) — S~{x,u, — ^) . Note that this special combina- 
tion (ID of twist-three amplitudes can be represented as n^^^^f^tws ^ j d^(j{-)(^^^ ^)6H^'"^{x, -0- 
It is free from singularities, since the integrand (5if*™^(x, — ^) = xH^^^ —^H^""^ is a continuous func- 
tion in the momentum fraction x. It remains an open problem if this cancellation of ln(— A^/ Q^) 
terms is ensured by general principles, or if it is only valid in leading order for the Wilson co- 
efficients and breaks down once radiative corrections are accounted for. Unfortunately, contrary 
to the case of general kinematics, in the WW approximation for DVCS we can not express T2 in 
terms of Ti. To test the WW approximation one has to use models for GPDs in order to study 
their influence on Ti and T2. 

The tree-level result for the angular coefficients up to twist-three accuracy in leading order of 
the coupling constant^ reads for 

• the interference term (p^) 

4 = -8(2 - y)xB^ |(2 - a;B)(l - y) - (1 - a;B)(2 - y)' (^1 - ^) | ReH, 

= -8{2 ~ 2y + y^)xBRen, = 8y{2 - y)xBlmn, (31) 



2 - xb 2 - xb 

the DVCS squared term (|16D 



eflf 



^The twist two contribution including gluonic transversity is known in next-to-leading order. 



11 



cj^'vcs _ 2(2 - 2y + y')xBnn* 



DVCS ?/)^B -p,„^eff-T/* „DVCS o t_ -lYefT-i 



„ ^ ReW^'H*, sr^^ = 8-^i^Im?i^"7i*, (32) 

2 — a^B 2 — xb 

where we defined the 'effective' twist-three function 

^es ^ (iT^'^ ^ ~ ^*"^) ^^^^ 

in such a way, that the remaining x-q dependence of the twist-three angular coefficients is the same 
as for twist-two. The coefficients cf and c^^*"^ contain, besides the gluonic twist-two transversity, 
also geometrical twist-four contributions. In leading order of perturbation theory the latter can 
be evaluated from handbag diagrams with additional transversal gluons and, so-called, cat-ear 
diagrams. As we observe the twist-three contribution induces a new cos(20) / sin(20) angular 
dependence of charge asymmetries as well as a new cos(0) / sin(0) angular dependence in the 
DVCS squared term with coefficients proportional to a 'universal' combination of GPDs given in 
Eq. (P3|). In contrast, the angular independent part of the interference term, i.e. cj, arises from a 
pure kinematical twist-three effect. 



4 Double and skewed parton distributions. 

Before we give examples for the shape of twist-three GPDs in the next section let us discuss the 
relation between skewed and double distributions (DDs). The double distributions were originally 
defined by § 

(P2rC+(K, -K)\Pi) = P+ dy dze-'^^+y-'^^+'F{y, z). (34) 

J-l J-l+\y\ 

By means of the so-called a-representation it has been shown that F{y, z) is a generalized function 
defined in the region |?/| -|- |z| < 1. This original definition, also introduced in [Q, leads to a 
violation of the polynomiality condition for GPDs 0. The reasoning goes as follows. The jth- 
moment of an GPD is given by the expectation value of local twist-two operators with spin j + 1, 
i.e. ^7?,^.^^ = '^pfjLi...^j^'^piVfjLi ■■■iVpj 'ip, which are completely symmetrized and traceless, 
and where V^=Vfj, — and = — igB^. Thus, they are polynomials of order j -|- 1 in 77: 

•'^^ k=0 
On the other hand 

/I r^^\y\ /*! dy I x — y\ 

dy dzS{y + r]z-x)F{y,z)= — S{y\x,r])F [ y, , (35) 
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with 



,1+r? ^ y ^ 1+ri 



+ eiv>x> -v)e f±2 > ^ > pi 



The jth-moment of H{x,t]) generates only a polynomial of order j in 77: / dxx^H{x,T]) = 
J dy J dz {y + rjzy F{y, z). As suggested in Ref. one can cure this problem by adding in 
the definition ( p^ of DDs an extra independent term concentrated in y = and proportional to 
A+: 

H{x,ri) = J dy J dz 5{y + rjz — x) {F{y, z) + ri5{y)D{z)} . (36) 

Since D{z) has the support \z\ < 1, it induces a term in the GPD that is entirely concentrated in 
the 'exclusive' region. Indeed this choice is unique and the D-term can be extracted from a given 
distribution by D{x) = limrj^oo H {rjx , 1]) . The first term in the bracket of Eq. (|36|) is understood 
as a representation of a subtracted GPD, namely, H{x,7]) — sign[rj)D{x/ri). Note that D{z) is an 
antisymmetric function in rj since H{x, —if) = H{x, r]). 

In the following we give an alternative solution of the problem stated above, in which the 
spectral function does not depend on the skewedness parameter 77. The new representation will 
enable us to derive an inverse transformation in a simple way, since it does not explicitly depend 
on 1]. It is instructive to start with local operators. The parameterization for the matrix elements 
of the completely symmetrized and traceless local vector operators sandwiched between spin-zero 
states reads: 

(^2|Xm...M. 1^1) =p,S, jP, . . . P,,H,^i,+i + A,P,, . . . P,,H,^i, + . . . + A, . . . A,^.i7,+i,o} . 

Next we introduce a generating function for the coefficients Hj^. It is easy to check that after 
contraction with the light-cone vector and summation over the local operators, i.e. 0+(fi;, —n) = 
UpTLl^K, -k), TZI{k, -k) = ^-^f^'^Mi • • • ^n'^lm-H' following definition of H{x, r]) 

1 



j dx x^ ^H(x,r])\r,=o, where < A; < j, 1 < j, 



is equivalent to that one in Eq. ([23| ) . Let us remark that the term proportional to Ap in the matrix 
elements of the twist-two light-ray operator 7?.^(k, — k) arises from combinatorics and is fixed by 



a WW type relation [|T5| 



(P2^7^J(/^,-/^)|Pl) = f^dxe-'^^^'' (^P,H{x,r]) + j\y W^ix , y) -^H {y , r])j , (37) 
where the kernel reads W2{x, y) = Of^{x, x — y). 
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Analogous to the definition of tlie GPD as a generating function of moments Hj{ri), we now 
introduce a double distribution f{y,z) with the moments Hj j^k = (^) I dy J dz y^~^z^ fiy^z) 
where < /c < j, 1 < j. Summing the series of moments up again we find 

H{x,ri) = J dy J dz x6{y + rjz — x)f{y, z). (38) 

so that the transformation (^) is modified in a minimal way by an additional factor of x. The ma- 
trix elements of twist-two light-ray operators are expressed in terms of the new double distribution 

by 

(P2^7^J(«:,-/^)|Pl) = JdyJ dzf{y,z) (|/P, + zA,) e-^'^^+^--^+^ (39) 



Now we are able to derive an inverse transformation for equation (P^. In the first step we 
project the Lorentz index p in Eq. ( p9|) onto the transverse plane and subsequently perform Fourier 
transformations w.r.t. A-|. and k. Employing the representation ( pT]) for the matrix elements in 
terms of GPDs, we immediately obtain the desired transformation: 

zf{y,z) = —-^-^ J dx J drj FV{y + rjz — x)'"^ J dx'W2{x,x')^H{x',ri). 

Here we used the following Fourier representation for |/«P+|: \a\ = J dX exp (iXa) PV^-, where 
we understand the PV prescription as PVp- = | { (A+'io)^ + (A-io)^ }• '^^^ ^^^^ steps are to perform 
the X integration and a partial integration w.r.t. t]. In the last step we dropped a surface term. 
We find 

fiy, = ^ / / {(y + '^^)~' -{y + ^v- xy^} H{x, ri). (40) 

Finally, we express the F— and Z)-functions of the representation (|36| ) in terms of the /- 
function. The D-term is easily extracted by taking the limit rj, x ^ oo with x/rj fixed: 

D{z) = z j dyf{y,z). (41) 

Consequently, the term J dy J dzx {6{y + rjz — x) — 6{riz — x)} f{y, z) corresponding to the first 
one on the r.h.s. of Eq. ( PB]) can be cast in the form ( PU]) by means of 

d f 

X {6{y + rjz — x) — 5{riz — x)} = y5{;y + rjz — x) + z— / dy'6{y' + rjz — x)W2{y', y)y. 

dz J 



It leads to 



F{y, z) = yf{y, z) - -^z j dy'W^iy, y')y'f{y', z). (42) 
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Thus, both F and D functions in (^) turn out to be different projections of the same function /. 
Note that the projection (|4^) together with the inversion formula (^) give us the function F{y, z) 
in terms of the GPD: 

F{y, z) = J dr]J dxPV [{y + zrj - x)^^} [H{x, t]) - sign{r]) D{x/r])] . (43) 

Here the D{x/t]) term appears from a surface term that arised from the convolution with the W2 
kernel [see r.h.s. of Eq. (42)]. 



5 Features of WW approximation. 

In this Section, which aims to serve illustrative purposes only, we will give qualitative estimates 
of the twist-three contributions. An analysis and numerical estimates for targets with spin will 
be given elsewhere. In the following we assume that the GPD factorizes in a partonic form 
factor F(A^) and a function depending only on the remaining three variables if(x, 77, Q^) = 
F{A'^)H{x, rjlQ"^). The non-trivial polynomiality requirement for the moments of H{x, rj) is obeyed 
by using the transformation formula discussed in the previous Section. To ensure the reduction 
of the GPD to the parton densities in the forward case, we assume that the double distribution 



can be modeled, following the proposal in Ref. |24], as a product of the quark densities with a 



profile function that has certain properties. In the case of the transformation (^) the ansatz for 
the quark DD is {y = 1 — y) 

F,{y,z) = My)7r{\ylz), vr(y, ^) = f,(y) = q^{y)e{y) - q,{~y)e{-y). (44) 

In the following we limit ourselves to the valence quark approximation, thus, we set g = and 
Q = Qvai which is, for the sake of simplicity, chosen to be determined by Regge and quark counting 
rules gvai(?/) = g« + Qd = ^V'^^^i^ - v)- 

In Fig. (a) the shape of the GPD that results from the transformation (|36|) with D = is 
shown (dashed). Taking a simple ansatzQ for the D-term, i.e. D{z) = 6{1 — \z\)2z{l — z"^), we 
see that such a term induces a complex shape (dash-dotted line). We also plotted the so-called 
forward parton distribution (FPD) model, where the GPD is equated to the quark density, i.e. 
7T{y,z) = S{z). As we already know, the WW approximation for the twist-three GPDs induces 
jumps at point |x| = C,- However, they cancel in the combination see Fig. |I| (b). 

In Fig. ^ (a,b) we show the DVCS amplitude Ti. as defined in Eq. ( ^2]) for the same twist- 
two models as introduced above. Here we observe that the DD-distribution model (|36|) without 
D-teim and the FPD model give almost the same DVCS amplitude. We also see that a D-term 

'*The sign and normalization are not relevant for our present discussion. 
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Figure 1: Twist-two GPDs H{x,ri) i] = \ are shown in (a) for the transformation (^) without 
(dashed) and with (dash-dotted) D term and the FPD model (dotted). In (b) we plot the corre- 
sponding WW approximation for 6H^'^^{x,ri) for ?7 = |. Changing r] to — ^ leaves (a) intact and 
reflects (b) w.r.t. the horizontal axis. 



can change the real part drastically, while it drops out in the imaginary part of the amplitude. 
It is worth mentioning that quite different shapes of the GPD do not change the simple convex 
(concave) form of the real (imaginary) part of the DVCS amplitude. 

Now let us address the WW approximation of the 'universal' effective combination ( ^31) that 
determines the normalization of the twist-three angular coefficients (|3T|), (0). We decompose it 
in WW and antiquark-gluon-quark pieces T-C'^iO = W^-"^"^ i^) + W^%0 . In the definition (^) 
we rewrite the derivatives^ and find 



n 



cflf-WW 



{H,{x,0-m-^,0}- (45) 



Since the integrand only has a logarithmic singularity at x = ^, the integration will result into 
a smooth ^ dependence, which, however, will be partly removed by the differentiation w.r.t. ^. 
As we can see in Fig. ^ (c,d) for the real part the normalization and shape of T-^^fr-ww ^jiffg^-g 
only slightly from that of Ti, while the differences in the normalization for the imaginary part are 
caused by the behaviour of the twist-two GPD at the point x = ^. In any case we observe no 
numerical enhancement of the twist-three contributions in the WW approximation. The curves 
suggest a direct connection of the WW approximation with the twist-two GPD in the small ^ 
region. For the FPD model we can even analytically perform the integration and the amplitude 
is entirely determined by the small x behaviour of the parton density. In the case q{x) = Ax^"" 
with a > for x — 0, we immediately obtain the imaginary and real parts: 



eff-WW 



1 



a /a 
+ if 2 



^Obviously, for any smooth function t{z), the following identity 



dx 52 — 



!(£/«) holds true. 
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Figure 2: The twist-two predictions (a,b) and the WW approximation (c,d) of 7^'^*^^^^, for the 
convolution of the hard scattering amphtude with the GPDs specified in Fig. (Q) are shown for 
the real (a,c) and imaginary part (b,d), respectively. 



ReH = -AirC" cot 
where C(a) = J^dxx~"' 



y) 



cff-WW 



n cot 



an 



aC{a) 



(46) 



ln(^)/(l + a;)-ln(i±^) /{1-x] 



It is important to note that a 



different ansatz can give quite different amplitudes in the small ^ region. This fact may be helpful 
in order to distinguish GPD models by comparison with small xb data from HERA experiments. 

Finally, we take a look beyond the WW approximation. Both the measurement of the polarized 
structure functions in deep inelastic scattering as well as lattice data suggest that in the forward 
limit the WW approximation is quantitatively valid. For the DVCS process we have to compare 
the size of the effective twist-three contribution ( |^ ) in the WW approximation with that one 
coming from the antiquark-gluon-quark correlation functions: 



i=u,d,. ^ 



dx-^\n 



2^ 



^ + x ^ — X — iO dx"^ 



Since the coefficient function contains now a double pole at x = ^, one has in general a numerical 
enhancement of S'j(x,^). Thus, the antiquark-gluon-quark contribution would have to be much 
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smaller than Hi{x,^) to favour the WW approximation. This can presumably be tested by ex- 
tracting different azimuthal angular moments from experimental data as discussed in the previous 
Sections. 



6 Conclusions. 

In the present contribution we have studied DVCS on a scalar target to twist-three accuracy. Let 
us summarize the lessons we have learnt from this exercise: 

• The interference term is dominated by cos0/ sin0 dependence away from the kinematical 
boundary. Q scaling of the Compton form factors can also be tested by a set of relations 
among the azimuthal angular coefficients. 

• Twist-three effects do not show up in the coefficients of the azimuthal angle dependence 
already present in the leading twist approximation but rather they induce new Fourier 
components. Thus, we expect at most (9(A^/Q^) and 0(M^/Q^) corrections to the leading 
twist angular dependence. 

• Twist-three functions contribute in a singularity free combination to the physical cross sec- 
tions and there is no violation of factorization. 

• GPDs are related to a single spectral density, or double distribution. 

• Rather different shapes of GPD models result in quite similar shapes for the predicted DVCS 
amplitudes manifesting low sensitivity to the former. However, the overall normalization 
differs significantly. 

A similar analysis for the DVCS cross section on the nucleon will be given elsewhere. 

We are grateful to V.M. Braun, M. Diehl, A. Freund, N.A. Kivel, M.V. Polyakov, A.V. 
Radyushkin and O.V. Teryaev for discussions on different aspects of this paper. 



Note added: Recently it has been pointed out in [g^] that the inversion of Eq. (^) for the 
function H{x,ri)/x is known as the inverse Radon transformation [26|. Our result (|^ ) differs 
from the standard formula by the first term in the integral on the r.h.s. of Eq. (^OD [and a surface 
term]. Obviously, our integral exist also in the case if H{x,ri)/x has a non-integrable singularity 
at X = 0. In the case when H{x,r])/x is integrable, the first term in the integrand can be safely 
neglected, being proportional to 6{y)6{z), and we obtain the inverse Radon transformation for 
the function H{x,ri)/x. Note that the PV prescription can be rewritten as PV J^^dz t{z)/z'^ = 
J^^dz{T{z)+T{-z)-2riO)}/z' ^ 
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